Fidelity and level correlations in the transition from regularity to chaos 
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Mean fidelity amplitude and parametric energy-energy correlations are calculated exactly for a 
regular system, which is subject to a chaotic random perturbation. It turns out that in this particular 
case under the average both quantities are identical. The result is compared with the susceptibility 
of chaotic systems against random perturbations. Regular systems are more susceptible against 
random perturbations than chaotic ones. 
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I. INTRODUCTION 

A. Peres wrote in 1984 a highly influential paper on the 
stability of quantum wave functions under random fluc- 
tuations of the Hamiltonian Qj. He introduced the over- 
lap between the wave-function propagated by a known 
Hamiltonian with the same wave function propagated by 
a slightly perturbed one as a measure for stability. This 
quantity is nowadays known as fidelity amplitude. Its 
modulus square is called fidelity or Loschmidt echoQ- 
Based on heuristic arguments, numerics and physical in- 
tuition, he concluded that a regular fluctuation, i. e. a 
fluctuations obeying the same superselection rules as the 
original Hamiltonian should have less effect on the sta- 
bility of the quantum state than a chaotic one, i. e. a 
perturbation with no additional symmetries. 

Later the same question was addressed amongst others 
[1, H| using ensemble theory [j| . The perturbation as well 
as the original Hamiltonian were chosen from a Gaussian 
random matrix ensemble (RME) . Average fidelity ampli- 
tude (f(t)) was calculated in second order perturbation 
theory. Exponentiating the result yielded 

ln</(t)) = -2„L (~ + j + C corr (r)) + G(r 2 ) , (1) 

where D is the mean level spacing of the original Hamil- 
tonian, t = t/tn is time measured in units of Heisenberg 
time tn = 2irh/D and T is the Breit-Wigner spreading 
width of an unperturbed eigenstate|6j. The first term 
is recognised as Fermi's golden rule (FGR). The second 
term is due to spectral fluctuations, i. e. fluctuations 
in the Hamiltonian which do not affect the eigenstates 
(called regular fluctuations by Peres). The last term is a 
correction term which accounts for the spectral correla- 
tions of the unperturbed system. 

Theperturbative result ([l} was confirmed experimen- 
tally |7H10| . Later it was completed and extended by ex- 
act calculations for different choices for the unperturbed 
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system and for the perturbation. Thereby unexpected 
features like fidelity revival at Heisenberg time [llj or fi- 
delity freeze [HHHl for purely off-diagonal perturbations 
were revealed. 

Fidelity is, at least when averaged over a large number 
of eigenstates, closely related to the parametric energy 
correlator, respectively to its Fourier transform, the so- 
called cross form-factor The latter measure the 
susceptibility of the spectrum against fluctuations of the 
Hamiltonian rather than the susceptibility of the wave- 
function. A relation between fidelity and these quantities 
is rather surprising but highly welcome for experimental 
purposes, since spectral measurements are much easier 
to perform than measurements of fidelity, which requires 
in principle knowledge of the entire wave-function. 

The above mentioned results were derived and are valid 
under the assumption that already the unperturbed sys- 
tem has chaotic dynamics. But the case originally con- 
sidered by Peres, where a regular system is perturbed 
by a chaotic admixture, is in quantum information de- 
vices more relevant. There the unperturbed dynamics 
is usually well controlled. It has been object of several 
theoretical and numerical d, HH-IH studies in the recent 
years. Therefore it comes as a surprise that for this case 
a detailed ensemble theoretical analysis akin to the ones 
performed in Refs.fTTl. IhH. fl5| is lacking. 

The present work aims at filling this gap. We consider 
the same situation as Peres did: A regular system is per- 
turbed by fully chaotic fluctuations. We present exact 
analytic results for the averaged fidelity amplitude and 
cross form-factor in the corresponding RME. It turns 
out that in average both quantities are identical! What 
is more the exact result is identical with the one obtained 
in exponentiated second order perturbation theory ([1}. 

The result is compared with the (known) results for 
originally chaotic systems, which are perturbed by ran- 
dom fluctuations, confirming that a regular system is 
more susceptible to random fluctuations than chaotic 
ones (Hill. 
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II. DEFINITIONS AND RESULTS 

Fidelity amplitude is defined by (h— 1) 

f(t) = ^\a im e- iH ^) . (2) 

Here Hq describes the regular system and the fluctuating 
Hamiltonian is given by 

H = H + \DV , (3) 

where V is a chaotic admixture. We average over the 
spectrum of Hq in an interval which is large enough to 
contain a large number N of (unperturbed) levels but 
small enough such that the mean level spacing D is con- 
stant. The strength of the perturbation is of order D. 
This means that the dimensionless perturbation strength 
A as well as a typical matrix element of V are of order 
one (see Eq. ((7|)). The parameter A is related to the 
Brcit-Wigner spreading width via T = 2irA 2 D. 

Following the work by Berry and Tabor [23| in a generic 
regular system, like for example a rectangular billiard, 
the energy levels are distributed in an intervall in the 
same way as independent random numbers. Assuming 
ergodicity the average over the energy intervall can be 
replaced by an ensemble average over N independent ran- 
dom numbers or likewise over an ensemble of N x N ma- 
trices with uncorrelated eigenvalues Em' , m — 1,...N 
(Poissonian spectrum). 

For definiteness we assume the distribution function 
w(Em^) of each eigenvalue of H to be a Gaussian with 
zero mean and variance N/2ir. In a region of order 
N around the origin the eigenvalues are uniformly dis- 
tributed with mean level spacing D = w(0) = N^ 1 ^ 2 up 
to corrections of order 1/N. This implies a weak form of 
translation invariance 

J dxw(x)f{x + y) = J dxw(x)f(x) +0(N- 1 ) (4) 

for any y = O(N ), which will be used frequently. 

We choose an incoherent superposition of all eigen- 
states in the interval as initial state. Equation be- 
comes 

/(*) - itre-^V*"* . (5) 

We define the cross form-factor as 

K(t) = l-tre- lHt tre lHot . (6) 

It is a purely spectral quantity and contains no infor- 
mation about the wave-function. This definition differs 
from to one of [l| by a singular contribution at t = 
and for N — > oo. 

The chaotic perturbation V is chosen from a Gaussian 
random matrix ensemble, defined by its second moments 

(VijVu)v = SuSjk + ( — - 1 J SikSji . (7) 
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FIG. 1: Plot of equation © for A = 0.1 (full line) and for 
a GUE perturbation. The curve is compared with fidelity 
amplitude of a chaotic original systems Ho, which is perturbed 
by a GUE with coupling strength A = 0.1. Here Hq is taken 
from the GUE (long dashed line, formula from [H]]), from the 
GOE (short dashed line, formula from [l^|) and from the GSE 
(dotted line, formula from [151]'). 

The Dyson index f3 labels the three classical ensembles 
[H]. The Gaussian unitary ensemble (GUE, p = 2, V 
Hermitean) is chosen when the perturbation apart from 
being chaotic breaks time reversal invariance (TRI) . The 
Gaussian orthogonal ensemble (GOE, /3 = 1, V real sym- 
metric) respectively the Gaussian symplectic ensemble 
(GSE, f3 = 4, V Hermitean selfdual) are chosen if the 
perturbation is chaotic but time reversal invariant. The 
GOE applies for integer spin and the GSE for half-integer 
spin. The variance of a typical matrix element of V is of 
order one. This means that the weak translation invari- 
ance (U) applies and the density of states of Ho remains 
unaffected by the perturbation. In the following angular 
brackets denote the above defined averages over both V 
and Hq. 

We now state our main results. First, in the limit 
N —> oo and for t > average fidelity amplitude and 
average cross form-factor are identical 

</(*)> = (K(t)) ■ (8) 

Second, for t > the exact result for the average fidelity 
amplitude is for all three ensembles 

(/(*)> = e-^(^) , r = ± . (9) 

EH 

This result coincides with exponentiated second order 
perturbation theory (JT]). For t — 0, (K(t)) differs from 
(/(e)) by a 5(£)-contribution. 

In Fig.Q]cquation © is plotted for a GUE pertubation 
(/3 = 2). The curve is compared with the corresponding 
ones for systems which are originally chaotic and which 
are perturbed by a GUE. We see that a regular system is 
less resilient against random perturbations than chaotic 
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ones. This is in accordance with former results [18j, l2l| 
in the literature. Among the chaotic systems the GSE 
is the most resilient against random fluctuations and the 
GOE the most susceptible. This indicates a monotonous 
dependence of fidelity decay on the spectral rigidity of 
the original system, as quantified by the parameter /3. 

III. CALCULATION OF FIDELITY 
AMPLITUDE AND CROSS FORM-FACTOR 

Despite their simplicity, the derivation of the main 
identities ([5]) and © requires a full fledged supersym- 
metric calculation. We sketch the main steps. 



Ay (J) and A^>(J) with respect to the entries of a super- 
matrix J. Finally the integrals over the vectors can be 
performed. 

After these steps both average fidelity amplitude and 
average cross form-factor can be written as (f(t)) = 
A / (J)Z(J,0|j=o and (K(t)) = A^(J)Z(J, i)| J=0 , where 
the generating function Z(J,£) has a supermatrix J as 
argument. It is given by the following supersymmetric 
matrix integral 

Z(J ' <} = Jrf d -^^- E ^i^E 2 ), (14) 

and 



A. Map onto a supersymmetric matrix model 

After Fourier transforming ([5]) and © fidelity ampli- 
tude and cross form-factor are expressed in terms of the 
resolvents G(E ± ) = (H - E ± ie)~ x and G (E ± ) = 
(H -E±ie)- 1 as 



/(*) 



1 

N 

k(t) = — 



dE x dE 2 j( El -E, 

dE x dE 2 i(Ei-E 2 
(2n) 2 



HtG q (E^)G(E+) (10) 
) *trG (^r)trG(S+) . 



Here and in the following we assume t > explicitely. We 
use the following fundamental property of the resolvent 



Gij(E) = 



1 d 



2dK, 



dot 



H -E + K 
H-E + K 



(11) 



K=0 



where K is a matrix containing source terms. Now both 
quantities can be expressed as derivatives with respect 
to source matrices K^ l > and K ^ of one and the same 
generating function 

J_ f dExdE 2 i( El -E 2 )t 
N J (2ir) 2 



iV J \AH j~ 

det(H Q -E l + JfW) det(H - E 2 + if (2) ) 
dct(H Q - Ei - KW) det(H -E+ - AT< 2 )) 



(12) 



as 



/(*) 
K(t) 



E 
E 

n , rn 



d 2 



(2) 
n i rn 



z 



k(V=o 



Jf(!)=0 

k( 2 )=o 



(13) 



After writing the determinants as Gaussian integrals over 
vectors with commuting (denominator) and anticommut- 
ing (numerator) entries the average over the perturba- 
tion (Z) v can be performed easily. The following stan- 
dard steps are explained in detail for instance in [25l — 
[27j . First a supermatrix is introduced via a Hubbard- 
Stratonovich transformation and then the differential op- 
erators in Eq. (1131) are expressed as differential operators 



Z(J, Ei,E 2 



d[a] exp 



r Str(cr — Ei) 2 )(15) 



2D 2 X 2 

(Sdet~ K (H <g> 1^ + 1 N ® E(J))) ff 



The square bracket (. . .)h denotes an average over the 
unperturbed Hamiltonian Hq. We have (n,p) = (1/2,2) 
for V £ GOE, (K,p) = (1, 1) for V € GUE and (n,p) = 
(1, 2) for V e GSE. Here E(J) is a Ap x 4p supermatrix 



E(J) = 



u~ 
E+l 2p 



(16) 



We use the standard definitions of a supertrace Strcr = 
tiAi — trA 2 and of a superdeterminant Sdetcr = 
dot A 2 1 dot (.4 1 — Ai J 4.7 1 A-2) of a block supermatrix a = 
' A x A x " 
A 2 A 2 

mmuting. The 2p x 2p supermatrix a has the form 



, where the entries of Ai and A 2 are anticom- 
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GOE 
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(17) 



The integration variables a n , n = 1 . . .4 are real com- 
muting variables and X n , n = 1 ... 2 are complex anti- 
commuting variables. The matrix integral extends over 
all independent entries of a. The 4p x 4p supermatrix 



J 11 
-J 21 



J 12 
— J 22 



(18) 



contains the source terms, where each of the 2p x 2p 
matrices J' J has the structure of the prototypes (jTTJ) . In 
a tedious but straightforward calculation the operators 
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A/(J) and A^-(J) can be worked out. They are given by 



2p / p 



A/ (J) (2p)2 E ( E dJZldJH 



nm mn 



A^(J) 




(2p) 2 



E 



(9 7 1 

n=p+l 
2p 



E 



^ 722 ^3 J 22 

^ ^ nn n— p+1 nn f 



(19) 



In the following we evaluate the matrix integral (|15p in 
the limit N — > oo, (£1 — E2)/D = r, and r finite. 



B. Evaluation of Eq. ([15]) 

The matrix integral over a is expressed in angle- 
eigenvalue coordinates a — > U^ 1 sll, d[cr] — > 
B(s)d[s]dfi(U), where U denotes the supergroup, which 
diagonalizes a. The diagonal matrix 

diag(s B i,isFi) , GUE, 
diag(s B i,SB2,«SFi,«SF2) , GOE, (20) 
diag(sBi,SB2,is F i,«SF2) , GSE, 

contains the Bosonic (s Bn ) and Fermionic ( isg p i) eigen- 
values. The Berezinian B(s) is given by [25t U 



B(s) 



1 



(s B l - ISFl) 2 
\S-B1 ~ gBaKfFl - SF2) 

(s B l - «SFl) 2 (sB2 ~ «SFl) 2 
IfFl " S F2 \S(s B1 ~ S B2 ) 

(is F i - SBl) 2 (iSF2 - S B l) 2 



GUE 
GOE 
GSE. 



(21) 



We proceed by observing that the average over Hq in 
the last line of equation (fl~5j) only depends on the 4p 
eigenvalues of E(J). These are only infinitesimally dif- 
ferent from the eigenvalues of E(0). The infinitesi- 
mal corrections are obtained by a perturbative expan- 
sion. The eigenvalues of S(0) are the eigenvalues of a 
and Ei which is 2p-fold degenerate. We lift this de- 
generacy by replacing E 2 t 2p with the auxiliary matrix 
E a = dwg(E% 1 ,...,E% p ,E% 1 ,...,E$ p ), such that E(0) 
is completely non-degenerate. Now the action of the op- 
erators Af (J) and A^-(J) can be mapped in a lengthy 
but straightforward calculation onto the action of first 
order differential operators in the eigenvalues of s and 



E a . We denote them by G f/K . They read 
G f 



1 ^ P( SB n,£g m )+P(sBn,-£g n 
(2p)2 



SBr? 



E 2 



1 £ V(ts Fn , fig J + V(is Fn , -E-J 



G K = 



{2pf 
1 



(2p) 2 



SBn - E 2 

]T V(s Bn ,-is Fn )V(E^ m ,E^ m ) ,(22) 



n,m 
P 



where we used the abbreviation T>(x, y) = d/dx — d/dy. 
The action of the differential operator with respect to the 
source matrix J on Z(J , E±, E 2 ) can be replaced by the 
action of G^ K as follows 



A f/R (3)Z(3,E 1 ,E 2 ] 



j=o 



(23) 



B(s)d[s]e-^ str(s - £l)2 G / /^Z (s,E a ) 



Here Zq(s, E a ) is the superdeterminant in the second line 
of equation (fTS"]) . The matrices s and E a are diagonal, 
thus the superdeterminant can be written as a ratio of 
ordinary determinants 



^o(s,E a ) 



Un=i d et(go ~ isFn) dct(g - E* n 
Tl P n=i<let(Ho-s Bn )dct(H -E*+) 



(24) 



Since the eigenvalues of Hq are uncorrelated the ensem- 
ble average over Ho of Zq is the iV-th power of a single 
integral 



Z (s,E a ) 



where R is given by 



dxw(x)R(x, s, E a ) 



iV 



(25) 



i?(x, S ,E a ) = (' n^i(^-^Fn)(a:--ggJ \ _ (26) 



Next we have to calculate the action of G^ K on Zq. We 
observe that at E a = E 2 t2 P the E 2 dependence of Z n 
drops out and Zq can be evaluated in the large A^-limit 



lim Z (s,E a ) 

N— >oo 



exp (-iimDNStrs) . (27) 



We need to calculate the action of G^' K on R. One finds 
in a straightforward calculation 



G s/R / dxw(x)R{x,s,V a ) 



, i?(s x) 

^£>(s B „,-js Fn ) / dx w(x) '-—^ , (28) 
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where 



R(s,x) = < 



(x — iSFl)(x — is F2 ) 

(x — s B1 )(x — s B2 ) 
x — ispx 



'Bl 



(a; — isFi)(x 



s m )(x 



GOE 
GUE 
GSE. 



(29) 



^B2> 



Equation (j28|) holds for both and for G K . From this 
follows immediately (f{t)) — (K(t)), which is our first 
main result ([5]). 

Collecting the former results and expressing everything 
in terms of the dimensionless energy difference r, we find 

(/(*)) = / TTT? e2mrT I d[s}B(s)e-^ Strs2 - l7I,iStIS 



(2k) 



P r 
^X>(s Bn ,-is Fn ) / dx w(x) 



R(s, x) 
x + r~ 



(30) 



We recall r = t/tn. Absolute convergence of the re- 
integration and of the s integration is guaranteed by the 
Gaussian weight functions. We are therefore allowed to 
interchange order of integration and perform the r inte- 
gration by the residue theorem. Moreover, we perform an 
integration by parts of the operator J2n=i ^( s Bn, — *SFn) 
using that this operator annihilates Strs and B(s). The 
result is 



(/(*)) = 2IL j d[s}B(s)e-^ strs2 -™ KStrs 

x Strs / dxw(x)e- 27rlXT R(s,x) .(31) 



The remaining integral is simple for the GUE but 
somewhat tricky for the GOE and for the GSE. 

In the GUE-case the ^-integration can be performed 
by the residue theorem and employing weak translation 
invariance. As a result the remaining integrals over sbi 
and over spi decouple. Both are Gaussian integrals. The 
final result is 



(/(*)> 



GUE 



(32) 



For the GOE the ^-integration is more complicated 
due to the square roots in R and the more complicated 
Berezinian B(s). The expression on the r.h.s. of (|3Tj) is a 
fourfold integral over the three eigenvalues of a and over 
x. In order to simplify this integral, we use identity (27) 
of Ref. [28| and thereby extract a GUE contribution from 
Eq. (1311) . Thus the average fidelity amplitude is given by 
the GUE result plus an extra term 

(/(t)>GOE = (/(t))GUB + (/(*)>add • (33) 

After using identity (28) of Ref. [28| and an integration 



by parts the extra term reads 

dsBldsB2dsFldsF2\sBl - sbsI 



</(*)> 



add 



2ir\ 2 J (sbi - «s F1 )(s B 2 - «s F1 ) 
1 



X(5(SF1 — SF2)exp [ -j-^jM 



Strs 2 - ^Str 



St 



rs 



^ T>(s B n, -iSFn) ~ 

dxw(x)e- 27riXT R(s,x) 



SBI - SB2 
— 2-jrixr ; 



■V(SB1, SB2) 



(34) 



The action of the differential operator in the third line of 
equation (|34|) on R(s, x) yields the crucial simplification 
of the integral 



</(*)) 



add 



2ir\ 2 



dsFids B ids B 2\sBi - S B 2\ 



rStrs" 



L Strs — 27rixr 



dxw(x)- 



(35) 



Li 2 



Now the (trivial) sfi integration decouples from the rest. 
Introducing coordinates v — sbi — sb2 and u = (sbi + 
sb2)/2 and using the weak translation invariance of w(x) 
it is seen that the ^-integration does not depend on u. 
Thus the w-integration decouples as well. Moreover the 
integrals over sfi and over u together yield again the 
GUE result. Thus we can write 



</(*)> 



add 



:</(*)) 



4ttD 
dxw{x) 



GUE 



dv\v\ 



e e 



\/x — v~ 3 \/x + v +3 



(36) 



The z-integration can performed employing weak trans- 
lation invariance of w(x) 

g — 2-Ki'XT ^ 

dxw(x) z — = — 2-KiD— 3q(2tttv) . 

yx — vr y/ x + v + 

(37) 

Here Jo (x) is the zeroth order Bessel function. After yet 
another integration by parts 

(/(t))add = (/(*)>GUE 

x ^-J (0)+ J dv\v\e~^J (2TTTv)j . (38) 

The remaining integral over v is a standard Gradsteyn 
integral With J (0) = 1 we find 



(/W)add = -(/(*)>GUe(] 



-2-k't^X' 



(39) 



and finally (/(t)) G OE = e^M^+i) . Again the re- 
sult coincides with exponentiated second order perturba- 
tion theory. Here it is even more surprising than for the 
GUE since the deceptively simple and compact result is 
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- at least in the way it was calculated here - the outcome 
of a complicated conspiration of terms. 

The GSE case is treated like the GOE case. The calcu- 
lation is simpler since no square root appears in R(s,x). 

The result is (/(t)) G SE = e" 4 ^^ ( . Thus we 
can write the mean fidelity amplitude in all three cases 
concisely as in equation (J9)). 

IV. CONCLUSION 

Using supersymmetry we calculated fidelity amplitude 
and cross form-factor in a random matrix model for a 
regular system with a chaotic perturbation. Surpris- 
ingly both quantities are identical under ensemble av- 
erage. The result is a simple exponential of two terms. 
One term decays quadratically in time and one term de- 
cays linearly. Exponentiated second order perturbation 
theory is exact, indicating that a more direct proof of our 
results is likely to exist. 

The fact that both quantities are identical can be rec- 
onciled with the general differential identity [HI, QjJ 

< /( '» = - W^)^' 8 ' (40> 

where only the subspace of the perturbation which is par- 
allel to the original Hamiltonian enters. To understand 
this we split the matrix space of perturbations V 3 V 
into a subspace V11 parallel and a subspace V_l perpendic- 
ular to the matrix space of the unperturbed Hamiltonian 



Ho 3 H via V G Vi trVHo = 0, ViJ € H and 
Vii = V — Vii . Likewise any perturbation can be written 
as V = V\\ + V± , where the parallel part shares the sym- 
metries of Hq- The model © might be generalized to 
H = Ho + DXu Vh + DX±V± where coupling strengths Ay 
and Aj^ might be different. 

Since a chaotic perturbation breaks any symmetry of 
the regular system in our case Vii consists in the trunca- 
tion of V to its diagonal part in the eigenbasis of H - It is 
intuitively clear and has been shown perturbatively [2fl ] 
that the diagonal part of the perturbation in the eigen- 
basis of Ho is responsible for the Gaussian decay. On 
the other hand the linear term (FGR) is due to the off- 
diagonal terms. This suggests to write our main results 
(0) and (J9) in the form 

(K(t)) = (/(*)) = e V J , (41) 

which obeys the differential relation (|40)) . Although we 
have proven equation (|41j) only for A|| = Aj_, we conjec- 
ture it to hold exactly for arbitrary coupling strength Am 
and Aj_. This conjecture is backed by a perturbative cal- 
culation akin to that leading to Eq. (fTJ). A rigorous proof 
remains a challenge for the future. 
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